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System parameters of the long-period polar V1309 Ori 
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Abstract. Based on high-resolution optical spectroscopy in the blue and the near infra-red spectral range, we 
derived velocity images (Doppler tomograms) of the mass-donating secondary star and the accretion stream of 
the long-period eclipsing polar V1309 Ori (RX J0515. 41+0104. 6). Combined with HST-spectroscopy of high time 
resolution and optical photometry we were able to derive the main system parameters and to determine the 
accretion geometry of the binary. The length of the eclipse of the white dwarf is At ec i = 2418 ± 60 s, the mass 
ratio Q = 1.37 — 1.63, and the orbital inclination i — 76.6° — 78.9°. The surface of the secondary star could be 
resolved in the Doppler image of Nal absorption lines, where it shows a marked depletion on the X-ray irradiated 
side. The accretion geometry proposed by us with a nearly aligned rotator, co-latitude <5 ~ 10° , tilted away from 
the ballistic stream, azimuth x — —35°, explains the shape of the emission- line Doppler tomograms and the shape 
of optical/UV eclipse light curves. 

Key words, accretion, stars: binaries: eclipsing, stars: individual: V1309 Ori, stars: magnetic fields, stars: novae, 
cataclysmic variables 
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1. Introduction 

Polars are cataclysmic variables (CVs), in which both 
stars rotate synchronized with the orbital period. The 
strong magnetic field of the white dwarf prevents the for- 
mation of an accretion disk. Instead it forces the accreted 
matter - when the magnetic pressure exceeds the hydro- 
dynamic pressure - to move along the magnetic field lines 
to one or both magnetic poles (see e.g. Cropper 1990). 

V1309 Ori (RX J0515.41+0104.6) was discovered in 
the ROSAT all-sky survey (Beuermann & Thomas 1993) 
and optically identified as magnetic cataclysmic variable 
by Garnavich etal. (1994) and Walter etal. (1995). It is 
one of the rare systems among the polars, which show a to- 
tal eclipse of the white dwarf. This gives restrictions which 
make it much more easy to derive the system parameters 
from observational data than for most of the other sys- 
tems. 

V1309 Ori has an orbital period of 7.98 h, by far the 
longest among all known polars. The spectral type of the 
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secondary star is between M0 and Ml (Shatter et al. 1995). 
A main sequence star of this type has a much smaller 
radius than the size of the Roche lobe at the given orbital 
period, i.e. V 1309 Ori is the only known polar with a 
significantly oversized or possibly evolved secondary. 

The magnetic field of the white dwarf could be infered 
from cyclotron harmonic emission to be ~ 61 MG (Shatter 
etal. 1995). There was some discussion in the literature, 
whether with the given field strength of the white dwarf 
and the long orbital period the system could rotate syn- 
chronously. Frank etal. (1995) showed that V1309 Ori 
can be understood in terms of the standard evolutionary 
model, if the system either was in an extended low state to 
bring the white dwarf into synchronism or that the mag- 
netic field of the secondary is strong enough to maintain 
synchronism. 

In the ROSAT light curves large short-term variations 
of the X-ray flux are seen attributed to the impact of in- 
dividual blobs of matter on the accretion region (Walter 
et al. 1995). 

Due to its exceptionally long period, V1309 Ori is 
highly suited for a detailed spectroscopic study involv- 
ing Doppler tomography. The long orbital period allows 
phase-resolved spectroscopy with high spectral resolution. 

In this paper we present high resolution spectroscopy 
of V1309 Ori, obtained with the two channel spectrograph 
TWIN at the 3.5-m telescope of the German-Spanish 
Astronomical Centre on the Calar Alto in Spain from 
1995/11/24 to 1995/11/28. Furthermore we present pho- 
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Table 1. The data sets of V1309 Ori used in this paper. (MLO: Mount Laguna Observatory, taken from Shaffer 
etal. (1995); DSAZ: Deutsch-Spanisches Astronomisches Zentrum, Calar Alto; AIP: Astrophysikalisches Institut 
Potsdam) 



Date 


Telescope 


Filter 


duration of observation 


time resolution 


1992/10/30 


MLO lm 


V photometry 


4:43 h 


30 s 


1992/12/23 


MLO lm 


V photometry 


3:20 h 


31 s 


1995/11/24 


DSAZ 3.5m 


spectroscopy 


4:56 h 


300 s 


1995/11/26 


DSAZ 3.5m 


spectroscopy 


3:05 h 


300 s 


1995/11/27 


DSAZ 3.5m 


spectroscopy 


7:44 h 


300 s 


1996/08/11 to 1996/10/26 


HST/FOS 


UV spectroscopy 


13 x 12 min 


0.8 s 


1999/11/04 


AIP 0.7m 


V photometry 


2:21 h 


30 s 


1999/11/16 


AIP 0.7m 


V photometry 


3:41 h 


30 s 



tometry taken with the 70-cm telescope at the AIP site. 
We also use V-band photometry published by Shaffer 
etal. (1995) and archived HST UV spectroscopy, which 
was recently used by Schmidt & Stockman (2001). 

From these data we derive the eclipse ephemeris, sys- 
tem parameters - inclination, mass ratio and the primary 
mass - and investigate the accretion geometry. 



Stockman (2001). This is time resolved spectroscopy in the 
UV with a time resolution of ~ 0.8 s. V1309 Ori was ob- 
served 13 times for approx. 12 min., including one ingress 
and four egresses of the eclipse. By phase-folding the data 
it is possible to create a composite eclipse light curve, cov- 
ering the phase interval from approx. -0.1 to 0.1. 

All observations used in this paper are listed in Table fy. 



2. Observations and reduction 

V1309 Ori was observed in 1995 at the 3.5-m telescope 
at Calar Alto with the double-beam spectrograph TWIN, 
covering wavelength ranges from approx. 4200 A to 5050 
A and from 7550 A to 8700 A. The spectral resolution of 
this data is -2.3 A (~ 84 kms" 1 at 8183 A) in the red 
channel and — 1.8 A (~ 115 kms" 1 at 4686 A) in the blue 
one. 

Flat fields, comparison lamp spectra and bias frames 
were taken regularly. 

For photometric calibration we took spectra of spec- 
trophotometric standard stars. Seeing and transparency 
were variable during our observations, the photometric 
calibration of the spectra is supposed to have an uncer- 
tainty of up to a factor 2. 

To avoid the long readout times of the CCDs, we put 
six - in slit direction shifted - single spectra onto the chips 
before reading out. This caused the night sky lines to be 
very bright in each spectrum, producing large errors at 
the corresponding wavelengths. 

By taking the data of all three nights, full phase cover- 
age is achieved with a phase resolution of 0.0125 (80 phase 
bins per orbital cycle). V1309 Ori could be observed for 
6.5 hours subsequently from Calar Alto. At the epoch of 
our observation, eclipses happened just before and just af- 
ter the phase of visibility, hence, eclipse phase is poorly 
covered by our spectroscopic observations. 

Photometric data in the V-band were obtained by ob- 
servations at the 70-cm telescope at the AIP in November 
1999. The observations with a time resolution of 30 s cover 
two eclipses. These data were calibrated with bias, fiat 
field and photometric standard star exposures. 

In addition we used an archived HST/FOS obser- 
vation, which was previously published by Schmidt & 



3. Results and analysis 

3.1. Determination of the ephemeris from photometric 
data (optical & UV) 

The knowledge of an accurate value of the binary period 
and of the inferior conjunction of the secondary star is 
a necessary prerequisite for any further modeling of the 
system. 

The long-term baseline for the determination of the 
period comes from the combination of our own photome- 
try with published photometry from Shafter etal. (1995) 
and an archival HST observation. 

The times of mid-eclipse of the new data were deter- 
mined using the same approach as Shafter etal. (1995). 
Straight lines were fitted to the ingress and egress parts of 
the light curves. The time of mid-eclipse was taken as the 
centre of the line drawn at half-ingress intensity between 
the two fitted lines. 



Table 2. The times of mid-eclipse used for the determi- 
nation of the ephemeris of V1309 Ori. 



Date 


Cycle 


time of mid-eclipse 


Error 


Filter 


1992/10/30 


-4250 


2448925.8348 


30 s 


V 


1992/12/23 


-4088 


2448979.7193 


30 s 


V 


1996/09/13 





2450339.43652 


15 s 


UV 


1999/11/04 


3449 


2451486.6147 


30 s 


V 


1999/11/16 


3485 


2451498.5890 


30 s 


V 
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Fig. 1. The V band eclipse light curves used for deter- 
mination of the ephcmeris and the composite UV contin- 
uum light curve. In addition, the UV line light curve is 
shown for comparison. In the UV continuum light curve 
the ingress and egress of the white dwarf are marked by 
dotted lines. 



By fitting a straight line to these times as a function of 
cycle number, the new photometric ephemeris was derived 
as 

HJD(T mid _ ecL ) = 2450339.4363(6)+£*0.33261194(8).(l) 



The numbers in brackets are the errors of the last digits. 

The eclipse data used is shown in Fig. |[ The phase 
for all data was calculated according to the ephemeris of 
Eq. (1). The UV continuum light curve was achieved by 
averaging the wavelength intervals 1319 — 1361 A, 1430 — 
1500 A, 1576-1612 A, and 1670-1698 A. For comparison 
purposes the UV line light curve, obtained by summing up 
the intervals 1381-1423 A, 1531-1569 A, and 1624-1659 
A and subtracting the interpolated continuum flux, is also 
shown. 

The deviations of the observed eclipse times from the 
calculated ones are shown in Fig. 




Cycle Number 



Fig. 2. The deviation of observed to calculated mid- 
eclipse times. Here the filled squares represent the 
data points used for our period determination (see 
Table ||). The circles denote the eclipse times from Shaffer 
etal. (1995), the empty squares the values published by 
Garnavich etal. (1994) and the triangles the values pub- 
lished by Buckley & Shaffer (1995). The cross marks the 
zero point of the spectroscopic ephemeris, derived from 
the velocity zero crossing of the Nal lines in our spec- 
troscopic data. The dotted line symbolizes the ephemeris 
from Shaffer etal. (1995) and the dashed line the one by 
Buckley & Shaffer (1995). 



As one can see, the orbital period from Eq. (1) is quite 
similar to the one of Shaffer etal. (1995). Our time base, 
however, is much larger than the one used there, making 
our period determination much more accurate. 

The ephemeris from Buckley & Shaffer (1995) (in- 
dicated by the dashed line in Fig. ||) predicts incorrect 
eclipse times for observations after 1995. 

Our spectroscopic observations allow us to determine 
a single value for the inferior conjunction from the radial 
velocity of photosphcric absorption lines. It is marked with 
a cross in Fig. [2] and is consistent with the determined 
period. 

The photometry presented by Shaffer etal. (1995) 
showed the eclipse length to be highly variable, which 
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means that not the white dwarf itself but the accretion 
stream contributes most to the observed flux. 



0.05 




-0.05 



2500 



Fig. 3. The UV trailed spectrum: The ingress and the 
egress of the white dwarf are clearly seen at the phases 
marked with small lines (-0.048, resp. 0.036). The phase 
is calculated from Eq. (1). 



The only data set, which allows us to identify the 
eclipse of the white dwarf, is the HST observation, which 
we show as a trailed spectrogram in Fig. | A UV contin- 
uum light curve based on these observations is shown in 
Fig. [j] (time resolution: 15 s). The data suggest that (1) 
the initial steep decline at phase -0.048 is caused by the 
eclipse of the white dwarf, and (2) that the eclipse egress 
of the white dwarf marks the end of the phase of totality 
(phase 0.036). 

The analysis by Schmidt & Stockman (2001) confirms 
this conclusion, but they claim the detection of a steep de- 
cline of length 6 ± 2 s, which they attribute to the eclipse 
of a small hot (> 150,000K) spot. Our analysis does not 
reproduce their results, we see the decline not only in the 
far UV but also in the near UV. We, therefore, attribute 
the decline to the eclipse of the whole white dwarf and 
consequently derive a much lower temperature for the 
structure which is eclipsed. Remaining line and contin- 
uum emission seen while the white dwarf is in eclipse, 
<fi = —0.048 . . . — 0.030, indicates that the accretion stream 
is not completely obscured during this phase interval. 

The ingress of the white dwarf eclipse lasts 45 ± 30 s 
(see Fig. |J), which is near the expected value of 71 s, if we 
calc ulate it with the system parameters derived in chapter 
3.3 and the mass-radius-relation for white dwarfs from 
Nauenberg (1972). 

An eclipse duration of the white dwarf of 2418 ± 60 s 
was determined (~ 0.084 phase units). Figs. || and |l| show, 




Fig. 4. The eclipse ingress in the UV continuum light 
curve with a time resolution of 15 s. The points of be- 
gin and end of the white dwarf's ingress are marked by 
dotted lines, while phase -0.0475 (dashed line) is equiva- 
lent to the center of the ingress. The phase is calculated 
from Eq. (1). 



that there is already some line emission ~ 0.005 phases (~ 
2 min) before the egress of the continuum light source. 

Eclipse ingress and egress of the white dwarf cannot be 
resolved in broad band optical photometry. This is due to 
the relative faintness of the white dwarf in that band - in 



Sect. 3.5 we estimate its contribution to only about 10% 
of the eclipse flux - and the rather poor time resolution 
of these observations of only 30 sec. 

As one can see in the UV continuum light curve in 
Fig. 0, the center of the eclipse of the white dwarf is not 
at phase zero, but approx. 0.006 phase units earlier. For 
conversion from photometric ephemeris to mid-eclipse of 
the white dwarf one has to subtract 172 (± 20) s from the 
time of mid-eclipse, so that the ephemeris of the white 
dwarf's superior conjunction is 

HJD(T wd ) = 24450339.4343(8) + E * 0.33261194(8). (2) 



3.2. High resolution spectroscopy 

A phase-binned trailed spectrum of the blue channel is 
reproduced in Fig. ||. It is clearly seen that the emission 
lines have different line shapes and intensities at different 
phases. It is possible to identify several components in 
the emission line. One can see that in the emission lines 
(e.g. He II A4686 and H/3 in Fig. ||) a narrow emission line 
component (NEL), nearly following a sine curve, and some 
more diffuse emission is existent. The NEL is supposed to 
originate from the irradiated hemisphere of the secondary 
star (see e.g. Schwope et al. 1997), while the other emission 
is likely to be originating in the accretion stream. 

In Fig. |^ we show the mean-orbital radial-velocity- 
corrected spectrum of V1309 Ori in the blue channel of 
the TWIN. For radial velocity correction we used a sine- 
fit to the NEL of He II A4686. 

In this spectrum several lines could be identified: the 
Balmer lines H/3 and H7, the lines of neutral (4387, 4471, 
4713, 4921 and 5015 A) and ionized (4541 and 4686 A) 
helium, the line of MgnA4481, the CIII/NIII-complex 
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Fig. 5. Trailed continuum-subtracted blue spectrum of V1309 Ori, created by phase-averaging the individual spectra 
of all three nights. 



6 

o 

DO 



2 40 




4200 



4400 



4600 4800 
Wavelength (A) 



5000 



Fig. 6. Mean-orbital spectrum of V1309 Ori in the blue channel, after correction to zero velocity (narrow emission 
line of HeiiA4686). 



at 4635-4650 A and two unidentified lines at ~4549 and 
~4584 A. These unidentified lines are interesting, because 
in the trailed spectrogram (Fig. ||) they show only the 
NEL component. 



Fig. |7| shows the phase-binned trailed spectrogram in 
the red channel. Some absorption and emission lines are 
evident while moving back and forth in wavelength. The 
area around the Nal doublet at 8183/94 A is magnified in 
Fig. | 



The spectrum shown in Fig. || was achieved by averag- 
ing the velocity-corrected (with a sine-fit to the Nal-lines) 
red spectra in the phase interval from 0.7 to 1.0. 

Several absorption lines (Nal AA8183/8194, KiA7664 
and KiA7699) and emission lines (HenA8236, OiA7775 
and the blend of Call and Paschen lines above 8400 A) 
could be identified. 

For further analysis it was necessary to determine the 
system velocity and to remove it from the radial velocity 
curves and trailed spectra. This was done by fitting two 
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Fig. 7. Trailed continuum-subtracted red spectrum of V1309 Ori, created by phase-averaging the individual spectra 
of all three nights. Emission features appear black, absorption features white. 
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Fig. 8. Mean spectrum of V1309 Ori in the red channel, after correction to zero velocity (Nal absorption lines) and 
averaging over the phase interval from 0.7 to 1.0. 



gaussians with fixed separation and same intensity to the 
Nal lines in each of 40 phase-binned spectra. 

By fitting a sine-curve to the centres of the gaussians 
as a function of the orbital phase, a system velocity of 
8(6) kms -1 was determined. With the baryocentric ve- 
locity correction of another 8kms _1 the (low) system ve- 
locity with respect to Earth of 7 = 16(6) kms^ 1 follows. 
The value of 8kms _1 was used to correct the observed 
spectrograms before tomographic analysis. 



3.3. The system parameters of V1309 Ori 

In this section we make use of the HST eclipse observation 
and radial velocity measurements of the secondary star in 
order to determine the orbital inclination i, the mass ratio 
Q = M w( j/M2, and the mass of the white dwarf M w d. 
As described in Sect. 3.1, the observed eclipse length 



of the white dwarf in the ultraviolet spectral range is 
At = 2418 ± 60 sec. For further analysis we assume, that 
the secondary star fills its Roche volume and that the 
white dwarf may be approximated as a point source. Then, 
according to Chanan et al. (1976), a purely geometrically 
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Fig. 9. Trailed spectra around HenA4686, around H/3 and around the Nal doublet at 8183/8194 A. In case of the 
latter the velocity is computed for the line at 8183 A. For clarity the data of Nal are shown twice. 



determined relation between i, Q, and Aip exists, which is 
shown in Fig. [Tc| . 




1.2 1.4 1.6 1.8 



MASS RATIO M,/M 2 

Fig. 10. For the measured eclipse duration of 2418 s the 
curve of possible combinations of inclination and mass ra- 
tio is shown (solid line). The dashed lines represent the 
possible combinations of mass ratio and inclination de- 
rived from the radial velocities of the NEL in our irradia- 
tion model for M 2 = 0.4 M Q (right) and 0.6 M Q (left). The 
dotted lines correspond to the uncertainties of the used 
parameters, e.g. the nearly vertical lines to the range for 
Ki and the other ones to the uncertainty of the eclipse 
duration. 



The parameters Q and i can further be constrained, 
if the orbital velocity K 2 of the secondary star is known. 




Fig. 11. The tomograms of Nal (left) and 
He II A4686 (right) with the outlines of the Roche lobes 
for Q = 1.42, M 2 = 0.6 M , i = 77.3° (dashed) and for 
Q = 1.56, M 2 = 0.4M©, i = 78.0° (solid). The contour 
lines are at 10% steps. 

Our spectroscopic data allow an indirect measurement of 
this quantity by observations of emission and absorption 
lines originating on the secondary star. According to their 
phasing and their photometric variation, the prominent 
narrow emission lines (NEL) must originate on the X-ray 
irradiated side of the secondary star. In principle, its ra- 
dial velocity variation gives clues to the orbital velocity 
of the secondary. The determination of the radial velocity 
amplitude is, due to severe line blending with components 
originating from the accretion stream, not done straight 
forwardedly. Recent experience gained from a study of an- 
other long-period polar, QQ Vul (Schwope et al. 2000), 
revealed Doppler tomography as appropriate tool to de- 
termine the NEL radial velocity. 

Doppler maps using Spruit's (1998) MEM code were 
computed for the principal emission lines H7, H/3, and 
He 11 A4686. The Doppler map of He 11 A4686 is reproduced 
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as contour plot in Fig. [Tl]. The radial velocity was deter- 
mined from Gaussian fits to an average profile along the 
i^-axis. All NEL- velocities agreed within 5kms _1 , sug- 
gesting a common location of emission. The velocity am- 
plitude used for further analysis is K 2 — 89kms . 

Also, the radial velocity variation of the Na-absorption 
lines can be used to constrain the orbital velocity of the 
secondary star. The trailed spectrogram of Fig. || imme- 
diately shows that the Na-lines have a much larger radial 
velocity amplitude than the narrow emission lines, that 
the radial velocity curve deviates from a sinusoid, and 
that the lines almost disappear at phase <f> — 0.5, when 
the irradiated side of that star is best visible. Altogether 
this means that the line formation and the radial veloc- 
ity curve is strongly influenced by irradiation from the 
accretion region on the white dwarf. This becomes obvi- 
ous from a Doppler map of the line doublet also shown in 
Fig. [ll]. The region near the inner Lagrangian point L\ is 
completely devoid of Na-line absorption. 

Radial velocities of Na-lines in individual spectra were 
measured by fitting a double Gaussian. The radial ve- 
locity amplitude of the lines were determined by fit- 
ting circular and elliptical orbit equations to the mea- 
sured phase-dependent velocities, we consistently arrived 
at K'i = 216 ± 7 km s" 1 . 

Using the measurements of the radial velocity ampli- 
tudes of the photocenters K 2 and K 2 of the irradiated 
and non-irradiated hemispheres of the secondary, the or- 
bital velocity can be calculated with an irradiation model. 
The details of the model were described just recently by 
Schwope et al. (2000) and need not be repeated here. It 
assumes that emission lines are formed at the irradiated 
part of the secondary's Roche surface and absorption lines 
only on the complementary region. For given period, in- 
clination, mass ratio and mass of the secondary factors for 
correction from centre of mass to centre of light result. 

The correction factor for the velocities of the photo- 
centers of the non-irradiated hemisphere to the center of 
mass, K2/K2, is only slightly dependent on the mass ra- 
tio. It ranges from 1.13 for Q — 2.5 to 1.18 for Q = 1.25. 
Hence, with the observed radial velocity amplitude of the 
Na-lines and the maximum correction factor the lowest 
orbital velocity of the secondary star is K 2 — 180 km s -1 . 
The resolution of the Doppler map of the Na-lines is not 
sufficiently high in order to decide whether the irradiated 
part of the secondary is at least partially visible in the 
Na-lines. We, therefore, allow also for a smaller correction 
factor K 2 /K 2 = 1.08. Our accepted range for the orbital 
velocity thus becomes K 2 = 180 — 200 km s -1 . 

In order to determine the mass of the secondary star 
in polars, one usually assumes that it is a main-sequence 
star. The size of the Roche-lobe is fixed for a given binary 
period. The mass of the secondary star can then be de- 
termined with the aid of a mass-radius relation for main- 
sequence stars. As pointed out by Garnavich et al. (1994) 
and Shaffer et al. (1995), the secondary is of spectral 
type M0-M1 and, hence, much less massive then a main- 
sequence star which fits in a Roche lobe at 8 hours orbital 



period (which would have about 1 M ) . The secondary 
star must have a mass as low as 0.6 M or less. Assuming 
M 2 = 0.6M Q , Kepler's third law gives Q = 1.16 - 1.38, 
and M wd = 0.70 - 0.83 M for the allowed range of K 2 . 
If we assume M 2 = O.4M , we get Q = 1.44 - 1.73, and 
M wd = 0.57 - 0.69 M . 

We also computed a grid of models for the NEL origi- 
nating on the irradiated hemisphere of the secondary star 
as a function of i, Q, and M 2 . The results are shown as 
iso- velocity lines for K' = 89 ± 5kms _1 for M 2 = 0.4 M Q 
and 0.6 M , respectively, in Fig. [HJ These lines cross the 
lines defined by the eclipse length almost orthogonally and 
give narrow constraints to the mass ratio and inclination. 
The accepted ranges are i = 76.6° — 78.9°, Q = 1.37—1.63, 
M wd = 0.82 - 0.58 M , and M 2 = 0.6 - 0.4 M . The out- 
lines of the Roche lobes for two combinations (the ones 
defined by the intersections of the dashed lines with the 
solid one in Fig. [h]) are plotted over the tomograms of 
sodium and HenA4686. The areas of origin of the spec- 
tral lines lie well within the Roche lobes. In principle, the 
size of the Roche lobe (i.e. Q) is determined by the ratio 
K 2 jK 2 = Vb/vf. A corresponding graph for Vb/vf vs. Q 
is shown in Schwope et al. (2000, Fig. 10). Our measured 
ratio i^Na/^NEL = 216/89 = 2.4 predicts a mass ratio 
Q = 1.63, at the edge of the range given above. We, there- 
fore, tend to accept a solution with a larger mass ratio and 
a smaller secondary mass. For further modelling the com- 
bination Q = 1.5, M wd = 0.70 M was used. 

3.4. Distance estimation 

The observed V-band flux of V1309 Ori during the eclipse 
is 17 m 45(10) (AIP, 1999/11/16) (in full agreement with 
Vcci = 17 nl 37 in Shaffer etal. 1995). We use YY Gem (A) 
as reference M0.5-star for an estimation of the distance. 
In Chabrier & Baraffe (1995) the absolute V-band bright- 
ness, the distance and the radius of this star are listed 
(M v = 8 m 99, R = 0.66 ± 0.02 R , D = 13.89 pc). 

Warner (1995) gives a formula for the Roche lobe 
equatorial radius i?i(eq) (in the orbital plane) as a func- 
tion of the mass ratio Q and the binary separtion a. For 
the most likely values in our parameter determination, 
Q = 1.5 and a = 1.48 x 10 9 m (for M 2 = 0.46), we cal- 
culate i?i(eq) = 5.15 x 10 8 m. This radius is 12% larger 
than the one of YY Gem (A) (4.59 x 10 8 m), which sup- 
ports the supposition of a larger radius of the secondary 
in V1309 Ori than of a main sequence star with the same 
mass. Both stars have the same spectral type, i.e. same 
surface brightness. Scaling the flux of YY Gem (A) to the 
larger area of V1309 Ori, one gets an absolute brightness 
of My = 8™74. The distance modulus then gives a reason- 
ably good estimate of the true distance of D = 550±75pc. 
The error includes the uncertainties of the mass ratio and 
of the secondary mass. 

Our distance determination has to be compared to 
those in the literature. Shaffer etal. (1995), by using the 
same method, arrive at the same distance (D > 500 pc) as 
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we do. They claimed their distance to be a lower limit due 
to a possible contamination by the white dwarf and the 
accretion stream during eclipse. The HST-spectroscopy 
shows, that the eclipse of stream and white dwarf are total 
so that the distance of 550 pc can be regarded as an exact 
value. 

Harrop-Allin etal. (1997) derive a much larger dis- 
tance of D > 1300 pc from K-band spectroscopy. This 
estimate is based on the non-detection of certain M-star 
features in a single spectrum which they interpret as due 
to the faintness of the secondary. However, their spectrum 
was taken during orbital phase 0.43 to 0.56, i.e. at phases 
where irradiation fundamentally changes the spectrum of 
the secondary as is evident from our phase-resolved Nal- 
spectroscopy. Hence, their estimate is not reliable. 

A value of D = 745 ± 70 pc, derived from the TiO- 
band strength, is given by Beuermann (1999). We cannot 
comment on the likely reason for the discrepant distance 
determinations, a possible cause might be non-solar abun- 
dances in the secondary's atmosphere. For further analysis 
we use D = 550 pc. 

3.5. Temperature of the white dwarf 

Based on our determinations of the mass M w d and the 
distance D we can estimate the temperature of the white 
dwarf from HST-spectroscopy. Due to low signal-to-noise 
in the individual HST-spectra it is difficult to extract a 
meaningful spectrum of the white dwarf from the data 
at ingress and egress phase. Instead we used the aver- 
age flux in the line- free spectral interval between 1430 A 
and 1500 A and compare it with model atmosphere spec- 
tra for pure hydrogen atmospheres kindly provided by 
B. Gansicke (Gottingen). The brightness steps at ingress 
and egress are not equal, that at eclipse egress is smaller. 
Whether the combined emission from the undisturbed 
white dwarf and a hotter accretion spot or an accretion 
flare right at eclipse ingress are responsible for this dif- 
ference cannot be answered with the present data due to 
the limited number of photons and due to the fact that 
this binary phase was covered only once with HST. In 
order to arrive at an upper limit for the white dwarfs 
temperature we take the flux difference at eclipse egress, 
which is -F1430-1500 < 4 x 10~ 15 erg cm -2 s^ 1 A -1 . Using 
-Rwd = 7 x 10 s cm for an 0.8 M Q white dwarf, which cor- 
responds to the white dwarf with the smallest surface 
within the derived mass range, and D = 550 pc we es- 
timate T wd < 20, 000 K. 

Schmidt & Stockman (2001) derived a lower limit of 
the white dwarf's temperature of 30,000 K. However, they 
made not clear how their value was derived. The reason 
for the differing results therefore cannot be identified. 

With our parameters the contribution of 
the white dwarf to the V-band is less than 
3 x 10 -17 erg cm -2 s _1 A -1 , which is less than 10% 
of the eclipse flux. 



3.6. The accretion geometry 

Clues to the accretion geometry can be derived from 
both photometric and spectroscopic data. In particular 
we make use of Doppler images of the main emission lines 
and the contact phases of certain features in emission line 
and continuum light curves. 
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Fig. 12. The tomograms of the emission lines: left (from 
top to bottom): H7 and Hen A8236, right: Hen A4686 and 
HeiA4471. Shown are also the Roche outlines for the 
above derived system parameters and the line for an as- 
sumed ballistic stream. 



In Fig. |l2j we show Doppler images of H7, He 11 A8236, 
HenA4686 and HeiA4471. The velocity resolution of 
the maps is of order 120 km s -1 within a radius of 
500 kms _1 around the origin. That the features in the to- 
mogram of H7 appear broader than the ones in the case of 
He 11 A4686 could be caused by blending with He 11 A4339. 
All lines, except HenA8236, clearly show the NEL from 
the irradiated hemisphere of the secondary. Beside the 
NEL, emission from the accretion stream can be recog- 
nized in the left and lower quadrants of the Doppler maps. 
The overlays show the size of the Roche lobe for our 



favourite parameter combination (Sect. 3.3) and the ex- 
pected location of the ballistic accretion stream. 

In all of our maps emission near the ballistic trajec- 
tory is visible. In He 11 A8236 we find emission only in the 
vicinity of the ballistic stream. It is possible that the devi- 
ation of this tomogram from that of He 11 A4686 is caused 
by the smaller signal-to-noise ratio. But in each of the 
other tomograms there is emission off the ballistic trajec- 
tory, which can only be due to emission from the magnetic 
stream. However, ballistic and magnetic stream do not ap- 
pear as separate structures in the map. It is, therefore, 
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not easy to determine a terminal velocity of the ballis- 
tic stream. We regard the spot of enhanced emission at 
( v x,Vy) — (—450, +30) km s -1 as the likely region where 
most of the matter couples onto magnetic field lines. 

Our tomograms are very similar to the ones of H/3 
and HenA4686 presented by Hoard (1999), where strong 
emission from the ballistic stream is seen. The only major 
difference is the claimed absence of the NEL in He II A4686, 
which is clearly present in our data. Possible explanations 
are the higher resolution of our data or a higher accretion 
state of the system during our observations, which caused 
a stronger irradiation of the secondary's surface. 

Although not identified as separate structure in the 
maps, some constraints on emission from the magnetic 
stream can be set. Firstly, emission is spread over a very 
extended region in Doppler coordinates indicating some 
kind of accretion curtain to be existent. Secondly, the bulk 
of emission is concentrated above a line v y = —200 km s _1 , 
which constrains the orientation of the velocity vectors 
along the accretion curtain. The maps differ in this respect 
clearly from those of e.g. HU Aqr, UZ For, and QQ Vul, 
which all show the magnetic stream to be bright in the 
lower left quadrant (Schwope 2001). These three rather 
well studied systems all have a 'standard' orientation of 
the magnetic field, i.e. the magnetic axis is inclined to- 
wards the ballistic stream. The discrepant appearance of 
the Doppler maps of V1309 Ori indicates a non-standard 
orientation of the magnetic axis. 

The HST-based emission line light curve shows stream 
emission up to phase —0.030, i.e. 0.018 phase units after 
the white dwarf is eclipsed by the secondary. Also, egress 
of the emission lines is observed prior to egress of the 
white dwarf by 0.004 phase units. Any model of the stream 
geometry has to account for this extended visibility. 

The optical B and V band light (sec Shafter et al. 1995) 
curves did not resolve eclipse ingress and egress of the 
white dwarf. Outside eclipse they are double humped with 
a primary maximum at binary phase 0.25. Around eclipse 
they show extended ingress and egress phases due to the 
eclipse of the accretion stream (or accretion streams). The 
start phase of totality varies between <f> = —0.024 . . . — 
0.037. The stellar photospheres of the white dwarf primary 
and the red secondary make small contributions to the 
BV light curves only, i.e. these bandpasses are completely 
dominated by the accretion stream. At the given inclina- 
tion and the extent of the ballistic stream, eclipse egress 
of the ballistic stream starts only at phase <f> — 0.068, 
when about 75% of the total light from the stream is vis- 
ible again after eclipse. We conclude, that the fraction 
of light originating from the magnetic stream/accretion 
curtain is about 75%, that from the ballistic stream only 
about 25%. The pronounced orbital variability finds an 
explanation in projection of an optically thick accretion 
stream/accretion curtain along the line of sight. Different 
height of the primary and secondary maxima in the broad- 
band light curves (e.g. m# = 15'. n 85 at phase <fi = 0.25 and 
m B = 16 m 30 at phase cj> = 0.65 on December 22nd, 1992) 
are explained by strong irradiation from the white dwarf. 



Obviously, the maximum projected area of the irradiated 
part of the stream/curtain is seen around phase 0.25. 

We developed a 3D model of the binary which predicts 
for a given set of binary parameters (i, Q, Ki, K2), for an 
assumed accretion geometry and for an assumed magnetic 
field strength the location of ballistic trajectories in spatial 
and Doppler coordinates, schematic Doppler tomograms, 
and light curves. For the generation of the light curves we 
assumed constant emissivity along the stream between L\ 
and the white dwarf. Hence, we cannot expect a proper fit 
to oberved data, but reproduce observed contact phases. 

With this model we interactively explored almost the 
whole parameter space within the limits on (i, Q, K%, K2) 
derived above. We assumed a dipolar magnetic field, it's 
orientation is described by the colatitude 5 with respect 
to the rotation axis and the azimuth (or longitude) x with 
respect to the line joining both stars. The quadrant be- 
tween x — 0° and 90° contains the ballistic stream. A 
satisfactory qualitative fit to the spectroscopic and photo- 
metric features described above was reached for 5 = 10°, 
X = —35°, i.e. with a nearly aligned dipole tilted away 
from the ballistic stream. Such a geometry in spatial and 
Doppler coordinates is reproduced in Fig. [l3[ 




Fig. 13. The accretion geometry, which reproduces the 
main features seen in the Doppler tomograms and the 
light curves. In the upper left corner the view perpendicu- 
lar to the orbital plane is shown, with the secondary star 
(right) and the accretion stream, falling onto the white 
dwarf (left). First it follows the ballistic trajectory, later 
it couples to the magnetic field and follows it to the mag- 
netic pole. The upper right panel shows the same sce- 
nario in velocity coordinates, plotted over the tomogram 
of He 11 A4686. The lower panels show the system at phases 
4> = —0.06 (left) and <j> = 0.068 with the correct inclina- 
tion angle. 
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Using these parameters we reach qualitative agreement 
between the observed contact phases in broad-band and 
emission line light curves, of the phasing of maxima in 
broad-band light curves, of the location of the stagnation 
region at the terminal end of the ballistic stream in spatial 
and Doppler coordinates and of the close vicinity of the 
magnetic and ballistic stream in Doppler coordinates. 

The predicted location of the accretion spot, which 
is offset from the magnetic axis, is at stellar colatitude 
17° and at azimuth of —16°. With this location of the 
spot we predict a maximum in the soft X-ray light curve 
(which is not yet observed) at binary phase <j> = 0.045. In 
the case of a pointlike accretion spot, a self-eclipse would 
be observed around phase = 0.55. If the accretion region 
is significantly extended on the surface of the white dwarf 
or in height, it will only be partially eclipsed. 

4. Discussion and conclusions 

The main results of our photometric and spectroscopic 
study with full phase coverage, paying particular attention 
to the eclipse properties is a re-determination of the sys- 
tem parameters of the binary. These are listed in Table §. 

Table 3. The system parameters of V1309 Ori, as derived 
in this work. 



Parameter 


Value 


P 


0.33261194(8) d 


^^ccl.wd 


2418 ± 60 s 


i 


76.6° - 78.9° 


Q = Mi/M 2 


1.37-1.63 


M 2 


0.6-0.4 M 




0.82 - 0.58 M 


K 2 


180 - 200 kms -1 


Kwd 


110 - 145 kms" 1 




< 20000 K 


D 


550 ± 75 pc 


S 


10° 


X 


-35° 



The orbital period derived here from broad band pho- 
tometry has much higher accuracy than previously, since 
a larger time base could be used by us compared to oth- 
ers. We found no evidence for an asynchronous rotation of 
the white dwarf. The time difference between the centre 
of eclipse of the white dwarf and of the accretion stream 
could be determined for one epoch by making use of an 
archival HST-observation to be 172 (± 20) s. 

A more accurate determination of the eclipse timings 
and of the eclipse duration of the accretion spot can prob- 
ably obtained from upcoming XMM-observations. Based 
on our geometrical model we predict no self-eclipse of the 
accretion region. We found no evidence of a second accre- 
tion region to be present but cannot certainly exclude its 
presence. 

Based on the measured velocities of the irradiated 
front- and non-irradiated back-side of the secondary star, 



we could narrow the possible range of masses of the sec- 
ondary to 0.4 — 0.6 M Q . This is the first dynamic mass de- 
termination of the secondary in V1309 Ori and confirms 
or even lowers previous mass estimates based on the spec- 
tral type. The secondary could be resolved in the Doppler 
image of Na absorption lines, where it shows a marked 
depletion on the irradiated side. 

Doppler maps of the emission lines are clearly different 
to those found in the literature for other eclipsing polars 
(HU Aqr, UZ For). We could explain these differences with 
an extreme orientation of the magnetic axis pointing away 
from the ballistic stream. 
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